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study was undertaken to examine temporal variation 








in blood- and environmental-lead levels in data observed for a 








sample of 249 children in Boston between 1979 and 1983 at the 








Brigham and Women's Hospital. The two primary objectives of this 


study were to: 








Determine th xtent to which blood-lead levels 
recorded in the study conducted at the Brigham and 
Women's Hospital exhibit seasonal variation. 








Determine if any existing seasonal trends in 
blood-lead levels are correlated with seasonal 
trends in environmental levels. 














For each child in the study, blood-lead and 








environmental-lead measurements were collected longitudinally 











over a period of two years. Levels of lead in air, dust, water, 








and soil were included in the environmental data. Nominally, 
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(blood or environmental lead) in six month increments. 
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was found to be more linearly related to blood-lead levels than 
was age itself. Consistent with other studies, blood-lead levels 
in children were found to increase with age. 

Air-, floor dust-, furniture dust-, and window sill 
dust-lead levels all exhibited highly significant seasonal 
variation. The estimated maximum-to-minimum ratios were 2.3 for 
air lead, -5 for floor dust lead, 1.4 for furniture dust lead, 
and 1.6 for window sill dust lead. Modeled lead levels for air, 





floor dust, 
peak window 


November. 


and furniture d 





sill dust 











ust all had peaks in July. 





ad 


levels wer 





Oddly, 


modeled to occur in 


Each of these responses were also significantly 





il 


related to the date of measurement, 








time. This is not unexpected 





the use of leaded gasoline. 


du 


with a decreas 


observed over 





to the concurrent reduction in 


The extent to which levels of 


lead in blood were 


correlated with levels of lead in the environment was also 





evaluated. As stated above, 





the seasonal 








in blood-lead levels was highly statistically significant. 


However, after adjusting for 








the 1 


inear 


l component of variation 


ffects of environmental 


measures, the (residual) blood-lead levels did not exhibit even 





marginally significant seasonal variation (at the 10 perce 





level). 





These results do not necessarily imply a causal 


relationship between seasonal 





vari 


ation 











nt 


in environmental-lead 


levels and seasonal rhythms in blood-lead levels. The fact that 


there were arguably parallel 
dust lead levels, doesn't imp 
influenced by the floor-dust 
floor dust-lead levels were t 








rhythms in blood- and, say, floor- 








ly th 





lead 





at the blood-lead levels are 


levels. In particular, 











o be 











retaining the same blood lead levels, an 
same statistical significance levels wou 
analysis approach. Thus, it would be im 


physiological model relating levels of 1 














those in blood, before proclaiming a cau 


Nonetheless in this data, 





early 1980's from a specific set of children in Boston, th 





abundant evidence supporting 
seasonal variations among blo 
furniture dust-lead levels. 


measures peak in July which ji 








month of June. In addition, 








environmental-—lead measures, 





same order of magnitude as the blood-lead ratio of 2.5. 


th 





which was collected in 





multiplied by two, while 


d models were refi 








ld be reported by 
portant to develop 


ead in the environ 





sal relationship. 











od-, 


Xistence and parallelism 


aur; 


floor dust-, and 





The three environmental-lead 


Ss very near 








the blood-lead pe 





the maximum-to-minimum ratios 


ranging from 1.4 to 2.3, 


a. 


if 


the 








bai A 





a 


the 


this 


ment to 


the 


ere was 


of the 


ak 





in the 
are of the 


Thus, 


based on the results of this study, it is quite plausible that 








seasonal variations in environmental-lead levels contribute to 


the blood-lead rhythms. 











iv 


TABLE OF CONTENTS 


















































Page 
1.0 INTRODUCTION 
2.0 DATA 
3.0 STATISTICAL MODELING APPROACH Sh ee ie Ye SE SD a a eo 3 
3.1 FORM OF SEASONAL VARIATION be) a ig es a we 4 


3.2 CORRELATIONAL DEPENDENCE AMONG 
REPEATED MEASURES ..... . 2. 0. ee ew ee ee 6 















































3.3 BLOOD LEAD ee ee ar a a NRE RE ee el is Ose tee tee te, Sas ee Se SS 9 
3.4 ENVIRONMENTAL LEAD SO wy te oe oe ae i SK we Le od 10 
3.5 BLOOD-LEAD LEVELS ADJUSTED FOR 

ENVIRONMENTAL LEAD LEVELS ............ 11 











4.0 RESULTS dB it a i eh a i A A ek Se AL 





DESCRIPTIVE STATISTICS : 
OUTLIER ANALYSIS god Bide Bove 
MODELING RESULTS FOR BLOOD LEAD ......... 
MODELING RESULTS FOR ENVIRONMENTAL LEAD ..... 2 
MODELING RESULTS FOR BLOOD LEAD AFTER 

ADJUSTING FOR ENVIRONMENTAL FACTORS ....... 22 





















































BEA BAD 
OB WNE 
































D0 DISCUSSBON 4% tice =H ao ek we oe he cake ome OAH Hs 26 














6.0 CONCLUSIONS AND RECOMMENDATIONS on ee Re wy ee Se ae we ee SOs 29 














REPERENGEH OS: Arf Ok ani ah A Ge 4 ee GPA Se A Se nid GP Binet UO 31 


LIST OF TABLES 


Table 1. Number of Observations Available for 
Each Response Across Time ........... 4. 2 
Table 2. Estimated Covariance Matrix for (Log) 








Blood Lead on an Individual Child ........ 8 

















Table 3. Estimated Correlation Matrix for (Log) 
Blood Lead on and Individual Child Wid Te Ro Ae Wo 8 

















Table 4. Geometric Means and Log Standard Deviations 
of Various Measures by Month eee & we & we we we 2 14 





Table 5. 


Multivariate Outliers 


16 





Table Mean Blood-Lead Concentration (pg/dl) by age 
(Controlling for Rate of Birth and Month of 
Measurement ) 

Table Results of Fitting Mixed ANOVA Model with Cyclic 
Seasonal Components to Blood Measures 

Table Results of Fitting Mixed ANOVA Model with Cyclic 
Seasonal Components to Environmental Measures 

Table Results of Fitting Mixed ANOVA Model with Cyclic 
Seasonal Components to Blood Measures Adjusting 
for Environments Lead Measures 

LIST OF FIGURES 

Figure 1. Unstructured covariance matrix 

Figure 2. Geometric average blood-lead levels with 95% 
confidence bounds by month and year 

Figure 3. Modeled seasonal variation and residual blood-lead 
levels after controlling for age and date of birth 
effects. (Bars represent 95% confidential bounds 
of blood-lead residuals.) 

Figure 4. Modeled blood-lead levels over time, showing effects 
of date, child's age, and seasonal variation J 

Figure 5. Modeled air-lead levels 

Figure 6. Modeled floor dust-lead levels 

Figure 7. Modeled furniture dust-lead levels 

Figure 8. Modeled window sill dust-lead levels 

Figure 9. Estimated seasonal component of blood and 
environmental lead levels, overlaid 

Figure 10. Blood-lead levels for five selected 








TABLE OF CONTENTS (Continued) 




















LIST OF TABLES (Continued) 

















































































































children 


born in January 


Vil 





21 


26 


les) 





23 


23 


24 


24 


25 


28 


Blood-lead levels for five selected 





Figure 11. 


children born in July 


viil 


28 


1.0 INTRODUCTION 





Several researchers have observed elevated levels of 
lead contamination and/or increased incidence of lead poisoning 


during summer months. Reasons for seasonal rhythms in blood-lead 





levels, if such a phenomenon is real, are not immediately 





apparent. The temporal variation may result from either altered 





2 





human physiology’’* or higher levels of lead exposure during the 


summer months. Determining the source of the temporal variation 








in blood-lead levels may enhance our understanding of the 





relationship between environmental-lead and its impact on body 








burden. 


There were two primary objectives of this study: 


e Determine th xtent to which blood-lead levels 
recorded in the study conducted at the Brigham and 
Women's Hospital exhibit seasonal variation. 














e Determine if any existing seasonal trends in 
blood-lead levels are correlated with seasonal 
trends in environmental levels. 

















This report examines temporal variation in blood- and 











environmental-lead levels in data observed on 250 children 


sampled in Boston between 1979 and 1983. 


2.0 DATA 

Umbilical cord blood samples were collected for 
11,837 births at the Brigham and Women's Hospital (formerly the 
Boston Hospital for Women) from April 1979 to April 1981. Of 

















these, 250 children were selected for an ongoing follow-up study 
involving environmental and psychological measurements. (One 
additional child had blood lead measured only at six months. 




















This data was used to estimate the average lead level at six 








months, but does not permit assessment of seasonal variation for 











this child.) The selection criteria for the follow-up study 








included cord blood levels in the highest, lowest, and middle 








deciles of the distribution of cord blood-lead levels, residence 


within 12 miles 


from the hospital, and likely 








to be available for 


two years of sampling. The resulting cohort differs from those 


studied in other research on lead exposure in 














that family incomes 





are relatively high, mothers were likely to be older, White, 
college educated, and working outside of the home. This analysis 
was based on the environmental and blood-lead data collected on 








each of these 250 children up to 24 months of age. 


For 





each child in the study, blood-lead and 





environmental-l 


The number of measuremen 
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analysis for blood lead and environment 





ad measurements were collected at various times. 
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in som 


Also, 





Cases 


they ar 


only referred to as lead 





information was not available regarding the 


proximity of these samples to children's activity areas. 


3.0 STATISTICAL MODELING APPROACH 


in this section. 


in 





soil dust. 


The approach to the statistical 


lead were blood, 


air, 


analyses is described 





The five media investigated for seasonal trends 





floor dust, furniture dust, and window 


Since there was a significant seasonal component in 


each response, it was also investigated whether a seasonal 





component remained present in the blood-lead levels after 





adjusting for the effects of lead in the environmental media. 





For the investigation of seasonal trends in each of 


the blood and environmental measures, each response was analyzed 





separately. When evaluating whether there was a significant 








seasonal component in blood-lead levels, after controlling for 





differences in environmental lead levels, the data were 


restricted to those sampling campaigns in which measures of lead 





in blood, floor dust, furniture dust, window sill dust, and air 


were obtained (see Table 1). This restricted the data to those 





collected at 6, 18, and 24 months. Excluding observations with 


incomplete data reduced the number of observations to 461 on 193 











children. The full data set contained 843 observations made on 


250 children. 








Fach blood and environmental lead measure was log 


transformed before analysis. There were two main reasons for 





this. First, these responses varied over one to four orders of 





magnitude. Second, after the log transformation all of the 


responses were better modeled by a normal distribution, which is 














an underlying assumption of the statistical analyses. 


Based on previous studies*’, factors suspected to 





influence children's blood-lead levels include the child's date 





of birth, age, and the time of year at which the sample is taken. 








These factors were included in the models fitted to the blood 








lead. For environmental lead, the dates of the measurement wer 








included instead of the child's date of birth and age to adjust 





for overall trends. 


3.1 FORM OF SEASONAL VARIATION 
The purpose of this analysis is to investigate the 


presence of a systematic seasonal component of variation in 








blood- and environmental-lead levels. It was not known at the 














outset whether such a cyclic component existed, let alone its 





functional form. However, it was assumed that a complete cycle 


for a seasonal variation, if present, would have a period of 





twelve months. 


The information in the data on date of sampling is 





limited to month and year of sampling. Therefore, the simplest 


and most general approach is to consider month as a class 





variable with twelve possible levels. This approach allows each 
month of the year to have its own mean after adjusting for other 


factors in the model: 





where x, denotes the row vector of covariates for each of the 
other factors in the model, $ denotes the column vector of 
parameters for the covariates, m(i) denotes the month number for 


the ith observation, ", denotes the deviation from the mean for 








12 
the mth month ( Ei", = 0} and ,; denotes random error. A 








limitation to this approach is the interpretation of the 12 


u 
. 
ml 


values of one is left with the burden of understanding 12 





different monthly averages. A second limitation is in the 





estimation of the variance of these monthly parameters. If a 





simpler model with fewer parameters underlies this cyclic 








variation, then estimates of its parameters would be more 


precise. 





Since for many of the media sampled, lead levels wer 





highest during the summer and lowest during the winter, a 
sinusoidal (Fourier) form of the model was investigated: 


Yi = x:$ + "cos ((m(i)-N)*2B/12) + ,4, 
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Figure 1. Unstructured covariance matrix. 


For example, F,, represents the covariance between 


taken at 12 months and 18 months on the same child, 





trolling for covariates. By definition F,, = F,, for 





,o,4s The diagonal terms, Fig, Foy» Fasp- Fao Cepreésent the 
for measures taken at each of the four increasing ages, 


ely, after correcting for other factors in the model. 





Without making any assumptions about the structure of 


jance matrix, Figure 1 represents the most general form 











This is referred to as unstructured. No relationships 








ferent covariances are assumed in the unstructured form. 


Below are two more structures which were considered: 





autoregressive and random child effect structures. 
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Autoregressive. The variance of blood-Pb 
concentration is assumed to be the same for 
children of all ages. The correlation 
between repeated measurements on the same 
child is assumed to decrease with time 
between measurements. 
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Random Child Effect. Variances are 
assumed equal for all measures. The 
correlations between repeated 
measures on the same child are 
assumed to be equal regardless of 
time between measurements. 
Observations on different children 
are assumed to be uncorrelated. 











ong these covariance structures is 





that although a more general model requires fewer assumptions 


about correlations, it requires th stimation of more 








parameters. To determine an appropriate covariance structure, 














the Akaike Information Criterion? (AIC) was used. The AIC 





provides a means of comparing models with the same fixed effects 











but different covariance structures. It is defined as 


where p(2) 


of parameters. 


IC = 2p(2)-2q, 





is the calculated log likelihood and q is the number 


This function adjusts the log-likelihood for the 





number of parameters in the model, including the terms for both 





fixed effects and covariance. The model having the largest AIC 


is considered to provide the best fit to the data. Due to the 








high AIC and the fact that variability was observed to decrease 








with the age of the child, the unstructured covariance was chosen 


as most appropriate. 





Table 2 displays the estimated unstructured 














covariance matrix for the model fit to the blood-lead data. 


TABLE 2. ESTIMATED COVARIANCE MATRIX FOR (LOG) 
BLOOD LEAD ON AN INDIVIDUAL CHILD 


Age 
(months) 








For illustrative purposes, Table 3 displays the correlation 





matrix associated with the covariance matrix in Table 2. 


TABLE 3. ESTIMATED CORRELATION MATRIX FOR (LOG) 
BLOOD LEAD ON AN INDIVIDUAL CHILD 








For consistency, no assumptions were made about the 
correlation structures of repeated environmental measures at a 
child's home. Therefore, unstructured covariance matrices were 


assumed for each media. 
The following sections describe the specific models 


fit to each of the responses. 


3.3 BLOOD LEAD 


Models of the following form were fit to blood-lead 




















concentrations: 
LB;, = $, + $, DOB; + $,a'/? + Siuay + > iar (2) 
where 
1 = child index, 
a = age of child in months, 
LB,, = the logarithm of the measured concentration 
of Pb in the blood (yg/dl) for the ith child 
at age a, 
$. = intercept, 
DOB, = date of birth of ith child, 
$, = limear effect of date of birth, 
$, = linear effect of age, 
Sin(ia) = seasonal effect for month in which child i 


was age a, and 


= random Error (ypas,- Pry, COtrrelated only af ay; 
i.e., measures are from same child). 





The model includes a date of birth effect to trace 
changes between different "birth cohorts" of children, and an age 
effect to reflect changes as a child grows regardless of his/her 


year of birth. To test whether there was confounding between the 











age and date of birth effects, the model was also fit without the 





date of birth effect. The estimates and significance levels 





obtained for the age effect in both cases were very similar. 





Therefore, it was concluded that the two factors were not 





confounded and both factors were included in our final model. 





The seasonal effect is described by the following 


Fourier model, 
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Sma) = 1008 ((m(ia)—-N,) *2B/12) + ",cos((m(ia)—-N,) *2B/6), 
where 
js = amplitude of annual cyclic variation, 


N, = phase of annual cyclic variation (time when 
peak occurs in months), 





; = amplitude of biannual cyclic variation, and 





N, = phase of biannual cyclic variation (time when 
first peak occurs). 








The two-phase cyclic model was selected as an objective 


compromise between a simple sinusoidal component with unknown 








phase, and a model including a different term for each month of 





the year. Certainly there was no reason to assume, a priori, 





that the shape of the seasonal component would fit a sine wave 











perfectly. A two-phase model was chosen by repeatedly adding 





Fourier terms with unspecified phase and period (12/k) months, 








k=1,2,3... until the relative reduction of error was 
insignificant. For blood lead this process was halted after 
adding the biannual cyclic term, which cycles twice per year, to 





the simple annual cyclic term. 


3.4 ENVIRONMENTAL LEAD 
Models of the following form were fit to levels of lead 


in air, floor dust, furniture dust, and window sill dust: 
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where the factors not defined above for the blood-lead model are 
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LA,, = the logarithm of the measured indoor air-lead 
concentration (pg/m?) in the ith child's home 
at age a, 


LFL, = the logarithm of the measured indoor 
floor dust-lead loading (pg) in the ith 
child's home at age a, 


LFU;, = the logarithm of the measured indoor 
furniture dust-lead loading (yg) in the 
ith child's home at age a, 





LWS;, = the logarithm of the measured window sill 
lead loading (yg) in the ith child's home at 
age a, 

eae = date when ith child was a months old, and 

$, = linear effect of date. 


The seasonal effect is described by the following 





Fourier model, for all four environmental media: 





aay = "cos ((m(ia)-N,) *2B/12) . 


m(ia) 


3.5 BLOOD-LEAD LEVELS ADJUSTED FOR 
ENVIRONMENTAL LEAD LEVELS 





A model was fit to the combined blood- and 











environmental-lead data to determine whether there were seasonal 





variations in blood-lead levels above and beyond those explained 








by changes (perhaps seasonal) in levels of lead in surrounding 





environmental media. The equation for this model is as follows: 
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where 
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ES 


LWA, 





Notice that the most general seasonal compon 


here, 


separately. 


the convention was to first 


seasonal 


after determining that the variation was significant 





controll 











variables 


formulation, 


the lo 


(ppb) 
ith ch 


the lo 


(ppb) 
child. 


as was the approach for fitting models 


(See equation 


but 


It is important to note 


in this model are subject 





environm 


esi 





not perm 


therefore, 





measurement 





ental lead measures). 





timates of these factors downward 


it an assessment of the magnitude of 


Crore: Thus, 


the measured 
soil outside 


garithm of 
of lead in 
ild, and 





the measured 
water at the 


garithm of 
of lead in 


1) 
m7 


ent, 














ing for other factors. 


that several of 





E to error (é€.9., Gach 


y 





These errors can pot 


(in magnitude). The 





it was not possible to adjust the estimates for 


it is reasonable to assume that o 


estimates of these effects are conservative. 


4.0 


This section presents the sta 


Sec 


describes the outlier analysis. 


Fit 





seasonal model to blood-lead level 





res 
in Section 4. 


fitting blood 


ting results. 


RESULTS 


tion 4.1 provides descriptive sta 


Section 





tistical analysis 











tistics. 


This is followed by the 











4.3 presents 


s. Corresponding model 





ults for the four environmental 


4. Finally, 


-lead level 














in Section 4.5 we discuss resul 





for environmental 


lead V 


lye 
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the predict 


these errors, 


Section 4.2 


concentration 
the home of the 


concentration 
home of the ith 


is used 


to each medium 


(1).) As mentioned in Section 3.1, 
t fit models with the most general 
to investigate simplification only 

t after 


COYr 





of the 


tentially bias 


data did 
and 


the 





ur 


results. 





mode] 


the results of fitting the 


l-fitting 


lead responses are discussed 





lts of 


s to the seasonal model after controlling 


4.1 DESCRIPTIVE STATISTICS 





For each of the responses measured there was evidence 


of periodicity. Table 4 displays estimates of the geometric mean 


levels (on original scale) of lead in blood, air, floor dust, 











furniture dust, and window sill dust by month. The number of 








blood samples collected is also listed by month. Approximate log 


standard errors of these estimates are provided for each media, 





along with the observed log standard deviations. The averages 


are often highest in the summer months (June, July, August) and 


lowest in the winter months (February, March). Four 


of the 





measures, blood, furniture dust, air, and window sill dust appear 


to have a relative minimum in September. For reasons discussed 














below, cord blood measures wer xcluded from the calculations 


for blood lead. 





Figure 2 displays observed geometric average blood-lead 


levels with 95 percent confidence bounds for each month of the 





study. These averages do not control for any covariates (such as 





age of child or date of birth). This figure reveals 





a slight 


cyclic variation, but does not show any sign of general change 








over time. It is not possible to distinguish between within- 





child effects (such as age) and between-child effects (such as 








date of birth) from this plot. The statistical modeling results 


presented next, allow this separation. It is shown 











within- and between-child effects actually counteract 


in this figure. 


4.2 OUTLIER ANALYSIS 


that the 





t each other 


Multivariate outlier analyses were performed to 


identify unusual data points. A Hotelling T* test was applied to 





identify potential outliers based on the distance be 


observation and the average of the remaining observat 








Cween an 





tions 


relative to the covariance matrix of the remaining observations. 
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The test was applied to the subset of the sampling campaigns in 





which each of the five measures: blood, air, floor dust, 


Hie) 


TABLE 4. GEOMETRIC MEANS AND LOG STANDARD DEVIATIONS 
OF VARIOUS MEASURES BY MONTH 


Number of 
Blood 
Measures 





pr 














Approximate 
Log Standard 
Error’ 


Log Standard 
Deviation” 











* 4 


The actual standard error of the mean log-transformed values varied due to 
sample size differences across months. These numbers represent the averag 
value of these log-standard errors. 














“ This represents the estimated within-month standard deviation of the log- 
transformed responses. 





furniture dust, and window sill dust were obtained. This test 


identified two observations as outliers at the 10 percent level; 














Table 5 displays the observations. The last column of the table 





provides the observed significance level of the Hotelling T* test 





for the observations. A Bonferroni-type critical value is used 














to compensate for the numerous simultaneous tests performed and 





to maintain a 10% overall significance. The second to last 





column is the appropriate threshold, based on the Bonferroni 





adjustment, to compare the observed significance level with. The 
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6th month floor dust measure on child 804391 was 600 pg, which 


was the largest value of floor dust lead when all 5 responses 
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100.00 


10.00 


1.00 


0.10 


Blood Pb Concentration (ug/dl) 


Dec 79 Jun 80 Dec 80 Jun 81 


Sampling Date 





Dec 81 Jun 82 Dec 82 Jun 83 


Figure 2. Geometric average blood-lead levels with 95% 
confidence bounds by month and year. 


were measured. The 18th month blood-lead measur 





was registered as 0 pg/dl (below the detection 











furniture and window sill dust-lead measures were h 


all observations with a 


fit with and without these two data points 








blood measure of 0 ypg/dl. 





effect on the conclusions. Since models fit 








to evalu 





tted to 





yielded the same conclusions as models fitted to th 





outliers removed, result 





including all of the dat 











ta. 
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on child 805874 
limit), but the 


ighest among 
Models were 
ate their 


all data 





e data with 


ts presented herein are based on analyses 


TABLE 5. MULTIVARIATE OUTLIERS 


Lead Levels 


Window 
Sill Dust 
Q 


Observed 
Significance 
Level 


Bonferroni 
10% Significance 
Cutoff 


Furniture 
Dust 
Q 


Floor 


Dust 
Q 


3.29 x 10° 
1.35 x 10% 


1.95 x 10% 
1.95 x 10% 


804391 
805874 


600 210 


12 


600 


18 160 





Lower order principal components were also used to 


visually inspect the data for outliers. The lower order 


principal components, by construction, are linear combinations of 





the factors with the smallest variance. Therefore, outstanding 





realizations of these principal components are often examined as 
There were no unusual observations noted in 


The 


tial outliers. 


of the fo 


poten 





a plot urth and fifth of five principal components. 











two outliers mentioned above were typical data points as measured 


by these principal components. 
.3 MODELING RESULTS FOR BLOOD LEAD 

Table 6 displays estimates of mean blood-lead 
concentration by age, adjusting for date of birth and month of 


measurement. 


th 


averag 


(Least-squares means are present 


ted which represent 





modeled mean for each age, 





observed level, 


which measurements were col] 





and averaging across 


lected.) 


holding date of birth at the 





the 12 months in 


The average blood-lead level 




















at 6 months was significantly less than those observed at the 
other ages, but cord blood ad levels were actually higher than 
those observed at 6, 12, 18, and 24 months (though not 


significantly). 








with the mothers' blood-lead level, and does not 





consistent with measurements taken at different 


blood-lead 1] 





vels wer 








xcluded from the subseg 
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Since cord blood lead may be more associated 


appear to be 


time points, cord 





uent analyses. 





The results displayed in Table 6 also suggest that the 








increase in blood lead has a nonlinear relationship with age. 











Since at best there are only four time points on each child, 
only one-parameter models were considered to fit this curvature. 
A model containing a term for the square root of age appeared to 


provide an adequate fit. 


TABLE 6. MEAN BLOOD-LEAD CONCENTRATION (pg/dl) BY AGE 
(CONTROLLING FOR RATE OF BIRTH AND MONTH OF 
MEASUREMENT) 


Age 
(months) Bound 











Table 7 displays the results of fitting the model to 





blood-lead levels. Age of child was found to be significant, 
date of birth was not; the magnitude of these effects is 
discussed below. As described earlier, both annual and biannual 
cyclic components of variation were used in the model to describe 


seasonal variation. Both components were sinusoidal. Phase is 





estimated in months. For interpretation, we assumed a phase of 





1.0 corresponds to January 15, phase 2.0 corresponds to February 








15, etc. Added together, the peak of the systematic seasonal 





component occurred in late June; lowest levels occurred in early 
March. The magnitude of this seasonal difference was 0.93 ona 


log scale. This corresponds to a multiplicative increase of 2.54 











in blood lead during late June over levels in early March on the 








same child. This number is calculated as the range of the 


seasonal component over a 12-month period. 
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TABLE 7. RESULTS OF FITTING MIXED ANOVA MODEL WITH CYCLIC 
SEASONAL COMPONENTS TO BLOOD MEASURES 


Std. Error 
rcept 0.792 0. 327, 


Ee ~0.010 0.4260 
ones 0.194 0.0005 
root) 








of Birth 











nnuat : 
mplitude) Veta 
1 N, oe 659 ig <0.0001° 
hase) (Jul. ean 
i nual 
mplitude) ei 0.264 0.094 
2 
i N, 5.999 0.381 0.0198 
(Jun 15, Dec. 15) (11 days) 


* Statistical significance of overall seasonal component. 
2 Significance of bi-annual cyclic component after controlling for annual 
cyclic component. 





















































Figure 3 displays the multiplicative factor 


corresponding to the modeled seasonal component of variance for 








each month of the year connected by a solid line. The bars 


overlaid on this plot display the mean and confidence bounds for 





the residuals of observed blood-lead levels after controlling for 





age and date of birth. There is clearly a seasonal component 











present in the residuals which parallels the estimated seasonal 








component. To facilitate quantification of the cyclic component, 





the values in this plot were scaled to force the minimum 








multiplier to be 1.00. Reference lines were placed at 1.00 and 








2.54, where the minimum and maximum occurred. The minimum 





appears in March, the maximum appears in June with a value close 


to that modeled for July. The values on this plot do not 








represent predicted lead levels, but rather the ratio of average 





lead levels for different months of the year to the level for the 











month with the lowest average levels, i.e. March (after 








controlling for age and date of birth). 
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Figure 4 illustrates the relative impact of each of the 


factors in the complete model fitted for blood lead. The figure 





covers a span of about three years from December 1979 through 
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Figure 3. 


Figure 4. 





10.0; 





Multiplicative Seasonal Variation 








0.1: 





Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov’ Dec 
Sampling Month 


Modeled seasonal variation and residual blood-lead 
levels after controlling for age and date of birth 
effects.(Bars represent 95% confidential bounds for 
blood-lead residuals.) 








— Born Dec 79 
-— Born Dec 80 
—--— Born Dec 81 











Blood Pb Concentration (ug/dl) 











1 T T T T T T T 


Dec 79 Jun 80 Dec 80 Jun 81 Dec 81 Jun 82 Dec 82 
Date 





Modeled blood-lead levels over time,showing effects 
of date, child's age, and seasonal variation. 
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December 1982. There are thr curves displayed in the plot. 





Each curve represents predicted blood-lead levels for a child 








born at a different time. The solid curve which begins earliest 








and represents the modeled blood-lead levels for a child born in 











December 1979. The most striking feature of this plot is the 





extensive cyclic variation about a generally increasing trend. 








This increase reflects the significant age effect. However, the 


cyclic variation outweighs the effect of age. Notice that for 











children born in December, although the estimated minimum 
seasonal variation occurs in March (Figure 3), when the age 


effect is added the relative minimum lead levels occur in 











February. 

Because data was collected for only two years on each 
child, the solid curve terminates after two years. The second 
and third lines represent modeled blood-lead concentrations for 





children born one and two years later (in December 1980 and 




















December 1981). There was a slight (and statistically 


insignificant) decrease in blood-lead levels with date of birth. 








This is reflected by the slightly lower starting points for 





children born later. 


4.4 MODELING RESULTS FOR ENVIRONMENTAL LEAD 
The final models fitted to the environmental media were 


described in Section 3.4. These models included a linear effect 





for date and a cyclic seasonal effect. An unstructured error 


variance matrix was assumed for the repeated measures at each 





child's home. Both the dat ffect and the seasonal effect wer 











statistically significant for each of the four environmental 





media investigated. 

Table 8 displays the estimated parameters for the fixed 
effects along with the standard errors and significance levels. 
Since the fitted model included only a single Fourier component 


for each media the phase listed equals the time, in months, of 
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the predicted maximum seasonal variation. For each medium, the 


significance of both parameters of the cyclic component is tested 
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TABLE 8 


Response 


Air-Lead 
(535 
Obs.) 


Floor 

Dust-—Lead 
(871 

Obs.) 


Furniture 

Dust-—Lead 
(872 

Obs.) 


Window 
Sill 
Dust-—Lead 
(863 
Obs.) 





simultaneo 


how close 


. RESULTS OF FITTING MIXED ANOVA MODEL WITH CYCLIC 
SEASONAL COMPONENTS TO ENVIRONMENTAL MEASURES 





Estimate 


pee 


Cyc 


1. 751 
-0.025 


ic Annual 0.408 


(Amplitude) 





Component 











6.815 
(Jul. 9) 


Cyclic Annual 
Component (Phase) 


ic Annual 
(Amplitude) 





1.903 


-0.008 


0.203 





Cyc 
Component 











6.896 
(Jul. 12) 


Cyclic Annual 
Component (Phase) 


ic Annual 
(Amplitude) 





1.359 


-0.012 


0.184 





Cyc 
Component 








7.349 
(Jul. 25) 


Cyclic Annual 
Component (Phase) 


clic Annual 
(Amplitude) 








2.823 


-0.012 


0.240 





Cy 
Component 











10.588 
(Nov. 3) 


ic Annual 
(Phase) 


Cyc 
Component 














Loge 


usly and was highly significant in all cases. 





[o) 


ile or 


aa 


oO 


ONS 


te 


[) 


eS 


et 


; ae] 





OV 2327. 
(7 days) 





Notice 


ther in time the maxima are predicted for air lead, 





floor dust 


lead, 








are predic 


ted 





and furniture dust lead. 


to occur in July. 
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Each of these peaks 


Figures 5 through 8 display model 


solid lines for air, 


floor dust, 





dust, respectively. 


geometric means by mo 


d lead levels as 








furniture du 


Overlaid on these plots 


nth with confidence boun 


and window sill 


th 


st, 





ar observed 


ds. Considering the 


numbers of observations represented, the models appear to fit 


well for so few parameters. 


in each response over 





time. 





observed for air-lead levels. 


Notice 





The most 


the sligh 





t decreasing trend 


drastic decrease was 


One obvious reason for this would 


be the coincident reduction in use of leaded gasoline in 


automobiles. 
Figure 9 di 
varia 


curve has been adjust 


direct 
seasonal components b 
flu 
the 
flu 





ctua 


environmental mea 




















ctuations of 


splays 


ed for 


etween 


sures, 


lead in floor dust, 


th 








trend effects. 








modeled season 


tion for the five responses overlaid for 


comparison of the phase and magnitude of 
the four modeled media. 


tions observed in blood were larger than 


al components of 
comparison. Each 


This figure allows 
the estimated 


The 








those observed in 


but were similar in phase to 








months after blood le 


ad. 


4.5 MODELING RESULTS FOR BLOOD LEAD AFTER 
ADJUSTING FOR ENVIRONMENTAL FACTORS 


furniture dust, 


and air. 


dow sill dust lead was predicted to reach its peak 4 to 5 


The extent to which levels of lead in blood were 


correlated with levels of lead in 


evaluated. Bach of t 





included in a model, 


and a class month effect to evaluat 


lead levels. 





the environmental mea 





h 


along with dat 


lead was significant before adjusting for 


Sures 


measured 





This model is described in det 


(with class mont 


it was not significant after adjusting for 


(p=.1148). 





The significant predictive fact 
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nvironment 


te of birt 





Whereas the seasonal component of variation 


th 


th 
t 





CO 


the environment was also 


al media were 


h and age effects, 


te seasonal rhythms in blood- 


in Section 3.5. 
in bl 
ffects of 

p=<.0001), 


hese effects 


ail 








ood 








linear 





effect, 





rs in this model were 


floor dust lead, and age of the child. The effect of soil lead 





on blood lead was marginally significant. Recall that only 
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Modeled air-lead levels. 


Figure 5. 
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Dec 79 Jun 80 Dec 80 Jun 81 Dec 81 Jun 82 Dec 82 Jun 83 


Sampling Date 


Modeled floor dust-lead levels. 


Figure 6. 
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Modeled furniture dust-lead levels. 


Figure 7. 
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Modeled window sill dust-lead levels. 


Figure 8. 
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Figure 9. Estimated seasonal component of blood and 


environmental lead levels, 


average levels of lead in soil 


child. 





Since the remaining environmental media 


observed as significant, a model was fit with age 


birth and the significant environmental factors. 


there remained a significant month 





smaller model, 





further investigation, 





lead ] 


overlaid. 


and water were available for each 


were not 
and date of 
Under the 


effect. Upon 


it was found that adjusting for air-lead 


vels reduces 








levels in addition to floor dust and soil 








the monthly component of variation 


marginally significant. 





final model fit. Thus, much, but not all, 


variation in these blood-lead level 











variations in floor dust lead, soil lead, 





eae 


from very significant to 
Table 9 displays the results of the 
of the seasonal 

s can be attributed to 


and air lead. 





One must realize that the fact that certain 








environmental factors were observed to be statistically 





significant in these models does not necessarily mean that there 





is a causal relationship. For example, just because floor dust- 





lead levels were observed as a highly significant effect in 














predicting blood-lead levels, one cannot conclude that seasonal 
increases in floor dust-lead levels cause corresponding increases 


in blood lead. 


TABLE 9. RESULTS OF FITTING MIXED ANOVA MODEL WITH CYCLIC 
SEASONAL COMPONENTS TO BLOOD MEASURES ADJUSTING 
FOR ENVIRONMENTS LEAD MEASURES 


Response Factor Estimate Significance 








[soir vent | fonts | aoe | o.0r7| 
[iced | | noe | aoe | 0409 | 
Month (12 level class variable) 


5.0 DISCUSSION 








The magnitude of the seasonal variation observed in 

















these data is substantial. The reader is reminded that for this 
study, blood- and environmental-lead levels were measured on 
these children every six months at best. To detect a cyclic 
component, one must observe levels which are systematically 














higher at one time during the year than at another. Since base 


blood-lead levels vary substantially across children, it is best 











to examine repeated measures on the same child. However, if the 


cyclic component were simply sinusoidal, the greatest within 
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child contrasts would be observed (in this study) for children 


born in 





of variation. For example, 


the month with the highest or lowest seasonal component 


if the seasonal variation was such 


that the lowest value occurred in January, and the highest value 


occurred in July, then children born in January or July would 








have the best chance of exhibiting major seasonal deviations if 








they were sampled at 0, 6, 12, 18, and 24 months of age. 


Figures 10 and 





11 illustrate this phenomenon for children born in 





January and July. Figure 10 displays the observed blood-lead 


level 











ls for five selected children born in January. Figure 11 


displays blood-lead levels for five other children born in July. 


The 











se figures illustrate the type of variation in blood-lead 





lev 











ls experienced by the children studied. Because measures 


were only taken six months apart, January and July were chosen 


because they portray the greatest seasonal contrasts. 





However, if a child was born in April or October the 











actual differences in the blood-lead levels at the times 
measured, due to seasonal variation, would be near zero. Plots 
of these levels would be more flat. These children would provide 











little added value in estimating the magnitude of seasonal 


var 


ayes 





jation. Moreover, sinc 


measures are taken six months apart 





is difficult to estimate 








the parameters of higher-order cyclic 


components. Thus, if feasible, monthly or quarterly measures 


wou 





ld provide much better information about seasonal variation 


from children not born in months where the maximum or minimum 


occurs. 


The reader is also reminded that the use of leaded 


gasoline was being phased out during the time this study was 


conducted. 





It is possible that leaded gasoline was the source of 





much of the lead in the environment - particularly in the air. 








Since more travelling is done in the summer than in the winter, 


changes in 


emissions from leaded gasoline may have been a 





significant contributor to the observed seasonal variations in 








environmental lead levels. 











Today, the use of leaded gasoline has 
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been virtually eliminated. Therefore, seasonal variations may be 








less pronounced today. 
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Figure 10. Blood-lead levels for five selected children 
born in January. 
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Figure 11. Blood-lead levels for five selected children 
born in July. 
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6.0 CONCLUSIONS AND RECOMMENDATIONS 
This report summarized an an 


in blood-lead levels for a longitudina 





alysis of seasonal rhythms 


1 study of blood- and 











environmental-lead levels on 250 child 
between 1979 and 1983. The following 
at: 





ren sampled in Boston 


conclusions were arrived 


@® There was evidence of significant seasonal variation 
in both blood- and environmental-lead levels observed 


in this study. 











@® Lead levels in blood, air, floor dust, and furniture 








in the late winter. 











dust were typically highest in the summer and lowest 


uch of the modeled cyclic seasonal variability in 





lained by adjusting for the 





ffects of environmental 


ead, specifically in floor 











M 
blood-lead levels was expl 
a ; 
d 


ust and air. 


The magnitude of the seasona 


blood lead was estimated to be 0.93 on 





corresponds to a multiplicative increa 





during late June over levels in early 








Thus for a child with blood-lead conce 











of 2.9 pg/dl, the predicted level in J 
7.4 po/dl. 





Tf such a seasonal component 


to exist in blood-lead levels today, i 





on the development of health-based sta 








warning levels. Specifically, it woul 











into account the month of the year in 





is at risk. The results of this study 


blood-lead threshold should be used in 





because a child with a marginal blood-] 





1 component of variation in 
a log scale. This 


se of 2.54 in blood-lead 





March for the same child. 


ntration measured in March 





une would be about 


of variation is confirmed 
t could have a major impact 
ndards and the setting of 
d suggest the need to take 
determining whether a child 
suggest that a lower 
February than in July, 


ead level in the winter is 





anticipated to have a much higher level 





in the summer. 
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Although the results of this study show evidence of a 











all lived in Boston in the early 1980s. 








seasonal variations can be very large in magnitude. 


before attempting 


levels considered 








to be a health 


risk, more current 


large cyclic component of variation in blood and environmental 
lead levels, the reader must recognize that the children observed 


The results suggest that 


However, 


investigations covering a broader population base over more 


varied geographic and socio-econ 


performed. Also, 


blood-lead variation, more frequ 


each child. Specifically, it would be of greater val 


to 





omic con 


better unders 














assessment of season 


children at more tim 


fewer time points. 





n 


al rhythms in blood lead to samp] 


e points than to sample more chil 
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ditions should be 


lue in 


le fewer 





ldren at 


to adjust health-based standards or blood-lead 


tand the nature of seasonal 
t measures should be taken on 


the 
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